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Work  accomplished  in  the  last  funding  period  is  summarized  below: 


(A)  Electronic  Properties  of  Semiconductor  Superlattices 

Semiconductor  superlattices  have  drawn  substantial  interest  in  the  last 
decade.  In  an  effort  to  understand  the  electronic  structures,  we  have  devel¬ 
oped  some  new  computational  methods,  including  the  muti-band  effective- 
mass  method,  the  Wannier-orbital  method7  and  the  bond-orbital  method15. 

We  have  used  these  methods  to  calculate  the  detailed  band  structures  of 
many  semiconductor  superlattices,  including  GaAs-AlGaAs,  GaSb-  AlSb,  and  , 

InGaAs-InP  systems.  The  muti-band  effective-mass  method  is  for  calculating  I 

quantum- well  confined  states  derived  from  bulk  states  near  the  band  extrema 
and  can  be  used  to  study  acceptors4’10  and  excitonic  states3,9,21  in  quan¬ 
tum  wells,  the  Wannier-  orbital  method  is  most  suitable  for  understanding 
the  inter-valley  scattering  in  conduction-band  states,  and  the  bond-orbital 
method  is  ideal  for  studying  the  mixing  of  states  in  superlattices  (or  coupled 
quantum  wells)  derived  from  near  zone-center  bulk  states. 

The  inter-valley  scattering  effect  is  important  in  short-period  superlat¬ 
tices.  For  example,  when  the  GaAs  layer  in  a  GaAs-AIAs  superlattice  is 
sufficiently  thin  (<  30A),  the  T-valley  derived  energy  levels  become  close  to 
the  X- valley  derived  levels,  causing  substantial  mixing  (T-X  mixing).  Such 
an  effect  has  been  observed  experimentally  by  applying  an  electric  field  to 
the  superlattice.  The  effect  may  have  useful  applications  in  electro-optical 
devices.  Mixing  of  T  and  L  states  can  also  exist  in  superlattices  involing 
materials  like  GaSb,  in  which  the  L-valley  is  the  lowest.  Using  the  Wannier- 
orbital  method,  we  have  examined  the  the  conduction  band  energies  of  (001) 
GaAs-AlrGai_r  As  and  AUGai-jAs-AlAs  superlattices  as  functions  of  layer 
thicknesses,  alloy  composition,  wave  vectors,  as  well  as  external  hydrostatic 
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pressure.7  We  found  that  there  can  be  substantial  mixing  between  the  T  and 
(001)  X-vaUey  states.  The  amount  of  T-X  mixing  is  extremely  sensitive  to 
the  layer  thicknesses.  Due  to  symmetry  considerations,  the  existence  of  such 
mixing  depends  critically  on  whether  the  barrier  material  (the  material  with 
higher  T  minimum  and  lower  X-minimum)  contains  an  even  or  odd  number 
of  atomic  layers  in  each  superlattice  unit  cell.  We  have  also  performed  de¬ 
tailed  theoretical  studies  of  the  L- valley  derived  states  in  (001)  GaSb-AlSb 
quantum  wells  and  superlattices.19  We  found  that  the  intervalley  mixing  be¬ 
tween  the  (1/2, 1/2, 1/2)  and  the  (1/2, 1/2, -1/2)  L-valleys  plays  an  important 
role  in  determining  the  electronic  structures  of  GaSb-AISb  quantum  wells. 

The  newly  developed  bond-orbital  method  was  used  to  calculate  the  va¬ 
lence  subband  structures  of  GaAs-Ga^^Al^As  quantum  wells  grown  in  dif¬ 
ferent  directions  (  [001],  [ill]  and  [110])20  and  Ini_xGazAs-InP  strained-layer 
quatum  wells.23  High  quality  superlattices  grown  along  the  [111]  and  [110] 
directions  are  available  only  very  recently,  Our  calculations  demonstrate  the 
differences  in  band  structures  for  quantum  wells  grown  along  different  direc¬ 
tions.  The  Inj-j-GaxAs-InP  strained-layer  quatum  well  has  important  appli¬ 
cations  in  quantum  well  lasers  for  optical  communication.  Our  theroretical 
studies  provide  some  guidelines  for  the  engineering  of  related  laser  structures. 

Most  superlattices  grown  to  date  consist  of  semiconductor  alloys  (e.g. 
GaxAh-xAs).  If  one  uses  an  alloy  as  the  well  material  in  the  superlattice, 
it  is  expected  that  the  disorder  effect  on  the  electronic  states  will  be  im¬ 
portant,  because  of  the  reduced  dimensionality.  We  have  developed  an  ap¬ 
proximation  technique  which  allows  us  to  efficiently  implement  the  coherent 
potential  approximation  (CPA)  for  treating  the  effects  of  disorder  in  alloy 
superlattices.14  Using  this  method  we  computed  the  density  of  states  of  su¬ 
perlattices  within  both  VCA  (virtual  crystal  approximation)  and  CPA.  Due 
to  the  two-dimensional  nature  of  the  bound  states  in  superlattices,  the  ef¬ 
fect  of  disorder  is  found  to  be  much  more  pronounced  for  alloy  superlattices 
than  for  bulk  alloys.  We  find  that  the  energy  shifts  due  to  disorder  effect  for 
typical  superlattice  states  are  around  10-100  meV. 

(B)  Optical  properties  of  semiconductor  quantum  wells  and  super¬ 
lattices 

We  have  developed  a  new  line-shape  theory  for  photoabsorption  in  super¬ 
lattices.  The  new  method  allows  one  to  accurately  and  efficiently  calculate 
the  photo-absorption  of  a  realistic  superlattice,  including  excitonic  effect, 
valence  band  mixing  and  inter-wellcoupling.  The  theory  predicts  realistic 
absorption  spectra  (without  adjustable  parameters)  for  GaAs-Gai_xAlxAs 
superlattices  which  are  in  exceelent  agreement  with  the  measured  spectra. 
It  is  the  predictive  value  that  makes  the  calculation  significant.  The  pred- 
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itions  include  some  important  physical  phenomena  such  as  the  line  shapes 
for  saddle-point  excitons  and  Fano  resonances  of  discrete  exciton  states  in¬ 
teracting  with  other  continuum  transitions.  We  have  applied  this  newly 
developed  line-shape  theory  to  calculate  the  magneto-absorption  spectra  of 
GaAs-Gai_j;AlxAs  superlattces.26  Our  theory  is  valid  for  any  strength  of 
magnetic  field;  thus,  we  can  follow  the  change  of  absorption  spectra  as  a 
function  of  the  magnetic  field.  We  found  interesting  change  of  line-shape  of 
the  exciton  resonances  related  to  band-structure  saddle-points  due  to  mag¬ 
netic  field.  Our  theoretical  spectra  are  in  very  good  agreement  with  recent 
magneto-absorption  measurements**1.  Our  calculations  are  useful  for  iden¬ 
tifying  the  origin  of  many  complicated  features  observed  in  the  magneto¬ 
absorption  spectra. 

Figure  1  shows  a  calculated  magneto- absorption  spectrum  and  the  cor¬ 
responding  experimental  spectrum.  As  seen  in  the  figure,  the  agreement 
between  theory  (with  no  adjustable  parameters)  and  experiment  is  quite 
good. 

(C)  Phonons  in  Semiconductor  Superlattices 

We  have  performed  detailed  theoretical  studies  of  phonon  modes  in  semi¬ 
conductor  superlattices  in  an  realistic  rigid-ion  model.  We  have  calculated 
the  dispersion  curves  for  phonons  propagating  in  any  directions,  taking  into 
account  the  long-range  Coulomb  interaction. 12,16,25  Our  theoretical  studies 
provide  the  understanding  of  two  key  features  of  optical  phonons  in  po¬ 
lar  superlattices,  namely  the  long  wave-length  anisotropy  and  the  macro¬ 
scopic  interface  modes.  The  theoretically  predicted  frequency  splitting  due 
to  anisotropy  is  in  good  agreement  with  the  available  data**2 

We  have  also  developed  a  phenomenonlogical  microscopic  theory  based 
on  the  rigid-ion  model  with  simplified  Coulomb  interaction.17  The  theory 
provides  some  integral  equations  which  allow  easy  calculation  of  the  disper¬ 
sion  curves  and  displacement  eigenvectors  for  long  wave-length  phonons  in 
superlattices.  The  results  obtained  by  the  simple  model  agree  well  with  those 
calculated  by  the  full-scale  rigid-ion  model.  This  simple  theory  allows  one  to 
gain  a  deeper  insight  into  the  intriguing  behavior  of  long  wave-length  optical 
phonons  and  the  interface  modes  in  polar  superlattices. 

(D)  Nonlinear  Optical  Properties  in  Superlattices 

Free  carriers  in  p-type  semiconductors  can  absorb  light  via  intervalence- 
band  transitions.  As  the  intensity  of  incident  radiation  increases,  the  optical 
pumping  rates  become  comparable  with  the  relaxation  rates,  and  the  carrier 
distribution  will  deviate  substantially  from  the  thermal  equilibrium  distribu- 
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tion.  The  absorption  will  be  saturated  at  sufficiently  high  intensity.  In  semi¬ 
conductor  superlattices,  inter-subband  transitions  have  nonlinear  optical  ef¬ 
fects  similar  to  those  from  bulk  intervalence-band  transitions.  Inter-subband 
transitions  can  happen  in  both  n-type  and  p-type  semiconductor  superlat¬ 
tices.  The  major  difference  between  the  two  is  that  for  n-type  superlattices, 
the  inter-subband  optical  matrix  elements  are  appreciable  only  for  incident 
radiation  polarized  along  the  growth  direction,  whereas  for  p-type  superlat¬ 
tices  the  optic  1  matrix  elements  are  appreciable  for  incident  light  polarized 
in  any  directions.  We  have  performed  calculations  of  the  saturation  behav¬ 
ior  of  the  inter-subband  transitions  in  p-type  semiconductor  superlattices6'24. 
The  nonparabolicity  of  the  valence  bands  due  to  band-mixing  effect  is  taken 
into  account.  Finite  carrier  lifetimes  due  to  the  carrier-phonon  scattering 
are  considered  within  the  deformation  potential  approximation.  Deviation 
of  hole  distribution  from  thermal  equilibrium  due  to  optical  pumping  is  cal¬ 
culated  by  solving  the  coupled  rate  equations.  We  find  that  the  linear  absorp¬ 
tion  spectrum  due  to  inter-subband  transitions  has  a  number  of  structures 
corresponding  to  transitions  from  the  first  few  subbands  (predominately  the 
first  heavy-hole  subband)  to  higher  valence  subbands.  We  have  also  stud¬ 
ied  the  absorption  spectra  of  a  probe  beam  in  the  presence  of  an  intense 
pump  beam.  The  pump  beam  can  cause  a  steady-state  non-equilibrium 
carrier  distribution,  thus  substantially  modifing  the  absorption  spectra  for 
the  probe  beam.  We  have  obtained  several  interesting  pump-and-probe  ab¬ 
sorption  spectra  and  they  may  have  possible  applications  for  light-by-light 
modulation. 

Change  of  index  of  refraction  due  to  the  applied  electric  field  (electo-optic 
effect)  has  an  important  application  in  providing  adjustable  phase  modula¬ 
tion  of  a  light  beam.  We  have  examined  theoretically  the  electro-optic  effect 
in  semiconductor  superlattices11.  Excitonic  effects  are  considered  for  wide- 
barrier  superlattices  and  are  found  to  substantially  enhance  the  electro-optic 
effect  near  the  onset  of  the  excitonic  absorption.  We  also  find  that  the  electo- 
optic  effect  in  small  band-gap  superlattices  (such  as  GalnAs-InP  systems) 
tends  to  be  larger  than  that  in  the  large  band-gap  (GaAs-AlGaAs)  superlat¬ 
tices. 

(E)  Indirect  Excitons  and  Electron-hole  Liquids  in  InAs  quantum 
wells 

Under  intense  photon  field,  the  two-component  electron-hole  plasma  is 
expected  to  exist  in  quantum  wells  in  various  phases:  electron-hole  plasma, 
electron-hole  liquid  (EHL),  exciton  gas,  biexciton  (exciton  molecule)  gas,  or 
Wigner  crystal.  In  order  to  observe  the  various  phases  of  the  electron-hole 
plasma,  the  lifetime  for  the  recombination  of  electron-hole  pairs  must  be 
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longer  than  the  relaxation  time  for  reaching  thermal  equilibrium.  This  could 
be  satisfied  in  semiconductor  quantum  wells  with  an  indirect  conduction 
band,  such  as  Ge-Sit-^Gex  systems.  We  have  proposed  a  rather  unique 
quantum  well  system  with  an  indirect  valence  band  which  may  also  be  a 
good  candidate  for  observing  the  electron-hole  liquid1 . 

With  the  presence  of  a  compressive  uniaxial  stress  in  a  small-gap  semicon¬ 
ductor  quantum  well  (e.g.  GalnAs-InP  systems),  the  first  light  hole  subband 
tends  to  move  upward  toward  the  first  heavy  hole  subband.  At  a  critical 
strain  the  energies  of  the  two  subbands  coincide  with  each  other,  and  the 
level  repulsion  between  them  can  cause  the  first  valence  subband  to  curve 
upward,  resulting  in  an  indirect  valence  band.  Thus,  the  lifetime  for  the 
electron-hole  recombination  become  long  enough  for  possible  realization  of 
the  electron-hole  liquid.  The  binding  energy  per  electron-hole  pair  in  the 
liquid  phase  will  be  enhanced  due  to  the  large  hole  density  of  states  at  the 
maximum.  We  have  studied  the  lifetime  and  energy  dispersion  of  the  indirect 
exciton  which  consists  of  an  electron  localized  at  the  zone  center  and  a  hole 
localized  at  the  maximum  of  the  indirect  valence  band13.  We  found  that  the 
lifetime  of  excitons  can  be  lengthened  by  three  to  five  orders  of  magnitude 
when  the  valence  band  goes  from  direct  to  indirect  under  the  uniaxial  stress. 
We  have  also  investigated  the  possible  formation  of  the  electron-hole  liquid 
in  In  As  quantum  wells  under  compressive  uniaxial  stress.  We  found  that  the 
electron-hole  liquid  phase  is  stable  against  the  exciton-gas  phase  within  the 
random-phase  approximation1. 

(F)  Qunatum  Monte  Carlo  Quantum  Monte  Carlo  Study  of  Mul- 
tiexciton  Complexes 

There  has  been  considerable  experimental  research  in  recent  years  R3  on 
bound  multiexcitonic  complexes,  systems  of  several  excitons  bound  to  shal¬ 
low  impurities  in  Si,  the  stability  of  which  is  enhanced  by  the  large  degeneracy 
of  low  lying  conduction  and  valence  band  states.  Recently,  evidence  has  also 
been  found  of  the  existence  of  free  complexes  of  three  and  four  excitons,  or 
polyexcitons  and  it  has  been  suggested  that  these  complexes  can  account  for 
a  recently  observed  unusual  condensed  plasma  photoluminescence  band  in 
SiR<4.  Theoretical  studies  of  bound  multiexciton  complexes  include  configu¬ 
ration  interaction  Rs  and  density  functional^6  approaches,  which  can  predict 
transition  energies  but  are  less  accurate  in  determining  ground  state  proper¬ 
ties  of  these  systems,  for  which  many-body  correlations  is  important. 

We  have  used  the  Variational  Monte  Carlo  method,  employing  a  Monte 
Carlo  technique  to  evaluate  variational  integrals  for  many-body  problems, 
to  calculate  the  binding  energies  and  ground  state  wavefunctions  of  two, 
three,  and  four  multiexciton  complexes,  in  the  spherical  effective  mass  ap- 
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proximation.  (  see  Figure  2  )  We  found  the  optima]  form  of  the  variational 
trial  wavefunction  to  be  a  linear  combination  of  Jastrow  many-body  wave- 
functions,  providing  a  description  of  a  system  intermediate  between  the  high 
density  electron-hole  liquid  and  a  complex  of  loosely  associated  excitons. 
The  results  for  smaller  complexes  have  been  compared  to  those  of  Green’s 
Function  Monte  Carlo  calculations  which  produce  exact  ground  state  ener¬ 
gies  within  the  limits  of  the  model  used;  the  variational  results  recover  about 
90  percent  of  the  binding  energy  of  these  systems.  Our  binding  energies  com¬ 
pare  well  with  experimental  results  for  Si1,2,  within  the  limits  of  our  model, 
although  central  cell  corrections  are  necessary  to  account  for  some  qualitative 
trends. 

FUTURE  WORK 

(i) We  will  study  the  atomic  geometries  and  electronic  structures  of  semicon¬ 
ductor  surfaces.  The  theoretical  studies  can  help  understanding  the  proper¬ 
ties  of  growth  surfaces  and  providing  guidelines  for  better  crystal  growth. 

(ii) We  will  study  the  electron-phonon  interaction  in  short-period  superlat¬ 
tices.  We  will  calculate  the  carrier  lifetimes  due  to  carrier-phonon  scatterings. 
We  will  consider  both  the  Frolich  mechanism  and  the  deformation-potential 
mechanism. 

(iii) We  will  study  the  transient  behavior  of  the  photoluminescence  excitation 
and  resonant  Raman  scattering  process  in  superlattices,  taking  into  account 
the  carrier  dynamics. 

(iv) We  will  combine  the  envelope-function  technique  with  the  bond-orbital 
method15  to  study  the  acceptor  states  in  quantum  dots,  wires,  and  superlat¬ 
tices. 

(v) We  will  continue  our  Quantum  Monte  Carlo  calculations  on  multiple- 
exciton  complexes,  including  the  central-cell  correction. 
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Figure  Captions 


Fig.  1.  Theoretical  absorption  spectra  (bottom)  and  experimental  pho¬ 
toluminescence  excitation  spectrum  taken  by  Deveaud  et.  al./n  (top)  for  a 
(oOA/  50A)  GaAs-Alo^Gao^eAs  superlattice.  Dashed:  heavy-hole  contribu¬ 
tion,  dotted:!ight-hole  contribution,  solid:  total.  All  theoretical  spectra  are 
broadened  by  1.5meV. 

Fig.  2  (a)Energy  of  N-exciton  complex  minus  energy  of  N  free  excitons,  in 
units  of  exciton  binding  energy  E(X),  as  function  of  the  electron-hole  mass 
ratio  <r.  Circles:  N  =  2;  X’s:  N  =  3;  Triangles:  N  =  4.  (b)  Energy  of 
N-exciton  complex  bound  to  a  donor  impurity,  minus  the  energy  of  a  donor 
and  N  free  excitons,  in  units  of  donor  binding  energy  E(D),  as  function  of 
the  electron-hole  mass  ratio  <r.  Circles:  N  =  1;  X’s:  N  =  2;  Triangles:  N  = 
3. 
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